suppressPackageStartupMessages(library(readxl))suppressPackageStartupMessages(library(tidyverse))suppressPackageStartupMessages(library(plotly))data =read_excel("Project_1_Data.xlsx", sheet ="pooled123")filteredData =select(data, PID, BSSQ_1:BSSQ_16,ASSQ_1:ASSQ_16,age, VRexperience, ssq_modelled)withVRexperience =filter(filteredData, VRexperience =='Yes')noVRexperience =filter(filteredData, VRexperience =='No')#calculating differences between baseline and active SSQ for each symptomfilteredData =mutate(filteredData, d_1 = ASSQ_1 - BSSQ_1)filteredData =mutate(filteredData, d_2 = ASSQ_2 - BSSQ_2)filteredData =mutate(filteredData, d_3 = ASSQ_3 - BSSQ_3)filteredData =mutate(filteredData, d_4 = ASSQ_4 - BSSQ_4)filteredData =mutate(filteredData, d_5 = ASSQ_5 - BSSQ_5)filteredData =mutate(filteredData, d_6 = ASSQ_6 - BSSQ_6)filteredData =mutate(filteredData, d_7 = ASSQ_7 - BSSQ_7)filteredData =mutate(filteredData, d_8 = ASSQ_8 - BSSQ_8)filteredData =mutate(filteredData, d_9 = ASSQ_9 - BSSQ_9)filteredData =mutate(filteredData, d_10 = ASSQ_10 - BSSQ_10)filteredData =mutate(filteredData, d_11 = ASSQ_11 - BSSQ_11)filteredData =mutate(filteredData, d_12 = ASSQ_12 - BSSQ_12)filteredData =mutate(filteredData, d_13 = ASSQ_13 - BSSQ_13)filteredData =mutate(filteredData, d_14 = ASSQ_14 - BSSQ_14)filteredData =mutate(filteredData, d_15 = ASSQ_15 - BSSQ_15)filteredData =mutate(filteredData, d_16 = ASSQ_16 - BSSQ_16)#reclasss VR experience as factor (was chr)filteredData$VRexperience =as.factor(filteredData$VRexperience)#we want to filter this data even further and split it into age groups#once in age groups, calculate the mean change for each of the age groups for each symptomfilteredData =mutate(filteredData, age_group =case_when( age >=16& age <=21~"16 to 21", age >=22& age <=29~"22 to 29", age >=30& age <=37~"30 to 37", age >=38& age <=45~"38 to 45", age >45~"above 45"))#convert the age groups into factorsfilteredData$age_group =as.factor(filteredData$age_group)str(filteredData)